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Abstract 

We present the complete next-to-leading order QCD corrections to the polarized hadropro¬ 
duction of heavy flavors which soon will be studied experimentally in polarized pp collisions 
at the BNL RHIC in order to constrain the polarized gluon density lS.g. It is demonstrated 
that the dependence on unphysical renormalization and factorization scales is strongly re¬ 
duced beyond the leading order. The sensitivity of the heavy quark spin asymmetry to lS.g 
is studied, including the limited detector acceptance for experimentally observable leptons 
from heavy quark decays at the BNL RHIC. As a further application of our results, gluino 
pair production in polarized pp collisions is briefly discussed. 


1 Introduction and Motivation 


Triggered by the measurement of the proton’s spin-dependent deep-inelastic structure function 
gi by the EMC [Q more than a decade ago, combined experimental and theoretical efforts 
have led to an improved understanding of the spin structure of the nucleon. In particular, we 
have gained some fairly precise information concerning the total quark spin contribution to the 
nucleon spin. The most prominent “unknown” is the elusive, yet unmeasured spin-dependent 
gluon density, Ag. Hence current and future experiments designed to further unravel the 
spin structure of the nucleon focus strongly on the issue of constraining Ag. In particular, 
information will soon be gathered for the hrst time at the BNL Relativistic Heavy-Ion Collider 

(RHIC) 0. 

The main thrust of the RHIC spin program is to hunt down Ag by measuring double 
spin asymmetries in longitudinally polarized pp collisions at high energies. RHIC is particularly 
suited for this task, since the gluon density is expected to participate dominantly in many 
different production processes. This is in contrast to deep-inelastic lepton-nucleon scattering 
where the gluon enters only as a small correction in the next-to-leading order (NLO) of QCD 
and indirectly via the renormalization group evolution of the parton densities. Along with the 
production of prompt photons and jets with high transverse momentum px, heavy flavor pair 
creation is one of the most promising candidates at RHIC to study Ag{x,pp) over a broad 
range of the momentum fraction x and scale pp- By virtue of the factorization theorem |^, 
helicity densities extracted from different processes should match. This important universality 
property, exploited in “global” fits of parton densities, can also be tested at RHIC. 

In the lowest order (LO) in the strong coupling ag, heavy flavor pair production in hadron- 
hadron collisions proceeds through two parton-parton subprocesses, 

qq QQ and gg QQ . (1) 

Gluon-gluon fusion is known to be the by far dominant mechanism for charm and bottom 
production in the unpolarized case in all experimentally relevant regions of phase space iii- 
This feature, true also in the polarized case unless Ag is exceedingly small, makes heavy 
quark production a particularly suited tool to study the gluon density. However, NLO QCD 
corrections to the LO subprocesses in Eq. (|I|) have to be included for a reliable description. First 
and foremost this is due to the strong dependence of the LO results on unphysical theoretical 
conventions such as the factorization scale, which reflects the amount of arbitrariness in the 
separation of short- and long-distance physics. In addition, the size of the NLO corrections 
turns out be quite sizable and not uniform in, e.g., the pp of one of the heavy quarks IB- 
The latter feature rules out the use of any approximations as estimates of the complete NLO 
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corrections. The computation of the NLO corrections are fairly involved since one has to keep 
track of the mass of the heavy quark, m, throughout the calculation. Massless approximations 
are bound to fail at small-to-medium values of pr where m ~ 0 {pt) and the cross section is 
large. So far only unpolarized NLO results had been available, see Refs. i0,i and @ for 
the differential and total cross sections, respectively. Polarized LO expressions can be found in 
Refs. P, P] , but the complete NLO results are presented for the hrst time in this work. 

Apart from calculational difficulties, a further complication arises when one tries to match 
theoretical parton-level results for heavy flavor production rates with experimental ones. Ex¬ 
periments can only observe the remains of heavy quark (meson) decays - usually leptons. In 
practice they also have to impose cuts on these particles to insure a proper c and b quark 
separation and to take care of the, usually limited and non uniform, detector acceptance. One 
thus has to hnd a practical way to incorporate hadronization, lepton-level cuts, and the de¬ 
tector acceptance in an analysis based on parton-level calculations, since they can distort spin 
asymmetries if polarized and unpolarized cross sections are affected differently. This may lead 
to incorrect conclusions about Ag. Heavy flavor decays usually have multi-body kinematics 
making it difficult if not impossible to trace back cuts to the parton-level analytically. Instead 
we propose to use “efficiencies”, to be dehned below, for bins in pt and pseudo-rapidity p of 
the heavy quark, to model its hadronization, decay, and crucial detector features. 


Heavy flavor production at RHIC is also interesting for reasons other than Ag. There is a 
well-known, longstanding discrepancy between data from the Tevatron collider for unpolarized 
open b production and theory |^. Recently, b rates in ep and 77 collisions were also found to be 
in excess of theoretical predictions • The fact that open c production seems to be fairly well 
described by theory makes these experimental Endings even more puzzling, since perturbative 
QCD should be more reliable for the much heavier b quark. This discrepancy for inclusive b 
production has revived speculations about new physics. For instance, the parameters of the 
minimal supersymmetric extension of the standard model can be tuned such that relatively light 
gluinos g exclusively decay into even lighter (s)bottom quarks b and 6 ’s: g ^ bb [|TI|. In this 


way the yield of 6 ’s is enhanced, and the Tevatron data can be reproduced while still complying 
with direct supersymmetry searches and precision measurements in e+e“ collisions [^. As 
will be briefly discussed below, our LO and NLO results for the dominant gluon-gluon fusion 
subprocess also contain the production of gluino pairs, gg —>• gg, after adjusting the color 
factors appropriately^ 0. Thus one can also estimate the spin-dependent hadroproduction 
of (light) gluinos at RHIC and at a conceivable polarized version of the CERN LHC in the 


distant future. Other explanations of the observed b excess make use of unintegrated 
gluon distributions and fc^-factorization, however, NLO calculations in this framework are still 
^The qq gg subprocess receives new contributions absent in qq QQ. 
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lacking. If results from the just upgraded Tevatron confirm the b excess, RHIC will play an 
important role in deciphering the underlying mechanism since it allows to study the energy- 
and spin-dependence of both c and b production. 


Finally, our results are also required for a fully consistent description of the polarized pho- 
toproduction of heavy quark pairs. Apart from the “direct” process 'yg QQ, where the 
NLO corrections have been calculated in Refs. [0, [T^, the (quasi-)real photon can also resolve 


into its hadronic content before the hard scattering takes place. The introduction of photonic 
parton densities is mandatory for the factorization of mass singularities of the direct process 
associated with collinear ■y ^ qq, q = u, d, s, splitting. Polarized “resolved” photon processes, 
Eq. (2), have been estimated [|^] to be small for fixed target experiments like COMPASS [R 


but can be significant at proposed future polarized lepton-hadron colliders like the EIC |17 


In this paper we focus on the outline and main results of our calculation. In Section 2 we 
give a brief description of the calculational steps and methods we have employed, in particular 
of the subtleties which arise due to the spin-dependence. In Section 3 we hrst discuss the differ¬ 
ent features of the NLO corrections to the polarized gluon-gluon fusion subprocess, always in 
comparison to the unpolarized case. Next we demonstrate the signihcantly reduced dependence 
on unphysical renormalization and factorization scales in NLO QCD for heavy flavor produc¬ 
tion at RHIC. We finish with studying the sensitivity of the heavy quark spin asymmetry to 
Ag at RHIC energies including realistic cuts on experimentally observable leptons from heavy 
quark decays. We conclude in Section 4. Technical details and lengthy analytical results are 
omitted throughout and will be presented elsewhere along with further phenomenological 
applications. 


2 Outline of the Technical Framework 

The 0{al) NLO QCD corrections to heavy flavor production comprise of three parts: the one- 
loop virtual corrections to the LO processes in Eq. (j^), the real “2 — > 3” corrections to Eq. (|I|) 
with an additional gluon in the final state, and a new production mechanism, gq{q) —> QQq{q), 
appearing for the first time at the NLO level. We choose the well-established framework of n- 
dimensional regularization, with n = 4-(-e, to tame the singularities of the loop- and 2 —> 3 phase 
space integrals. Ultraviolet singularities show up only in the virtual corrections and are removed 
by on-shell mass and coupling constant renormalization at a scale g,R. The latter is performed 
in a variant of the MS scheme which is usually adopted for heavy flavor production i,i0- 
This prescription is characterized by the decoupling of heavy quark loop contributions to the 
gluon self energy and leads to a fixed flavor scheme with nif = nf — 1 light flavors active in the 


3 





















running of and in the scale evolution of the parton densities. Infrared (IR) divergencies 
of the virtual diagrams are cancelled by the soft poles of the 2 —> 3 contributions. This includes 
also double, 1/^^, pole terms which show up when IR and mass/collinear (M) singularities 
coincide. The left over 1/e M singularities are then absorbed into the bare parton densities by 
the standard factorization procedure in the MS scheme. Other renormalization/factorization 
schemes can be obtained easily by additional finite scheme transformations. 


The required squared matrix elements |Mp for both unpolarized and longitudinally polar¬ 
ized processes are obtained simultaneously hj calculating them for arbitrary helicities Ai ^2 =“±” 
of the incoming quarks or gluons, i.e.. 


|M|^(Ai,A 2 ) = |M| +A 1 A 2 AIM/ 


( 2 ) 


using the standard helicity projection operators for bosons and fermions (see, e.g.. Ref. 

_2 

Results obtained for the unpolarized \M\ can be compared to the literature isia which 
serves as an important consistency check for the correctness of our new helicity dependent 
results A|Mp. To facilitate this comparison we closely follow the calculational steps and 
methods adopted in 1,0. It should be noted that, contrary to the unpolarized case [0], the 
processes qq QQg and gq QQq are not related by crossing for polarized initial states and 
have to be calculated separately. We should also recall here the dehnition of the spin-dependent 
parton densities, 

Af{x,fiF) = - f^{x,fiF) (3) 


where /+ (f^) denotes the probability to find a parton f = q, q, g at a. scale fip with momentum 
fraction x and helicity -1- (—) in a proton with helicity -|- [the unpolarized parton densities 
f{x,fiF) are obtained by taking the sum in (^]. In the following the compact notation 0 
denotes both an unpolarized quantity 0 and its longitudinally polarized analogue A0. 


The virtual (V) cross section for the qq and gg initial states is obtained up to 0{al) only from 
the interference between the virtual and Born amplitudes. Loop momenta in the numerator are 
dealt with by applying an adapted version of the Passarino-Veltman reduction program to scalar 
integrals which properly accounts for all possible n-dimensionally regulated divergencies 
in QCD. The required scalar integrals can be found in Qj, however, we have checked them 
by standard Feynman parameterization techniques (see also 0). The renormalized, color- 
averaged results can be decomposed according to their color structure as 
Cf 


\^qq\v — 9 


AN, 


c 


\Mm\l = 


2CfNqed + CaNqk + Nqf 
1 


2(iV2 - 1) L 


{ 2 Cf) Uqed + CaUqq + Ukq + 2 CfUrf + CaUql 


(4) 

(5) 


where g"^ = Anas, Ee = 1/(1 + ^/2), and AF^e = 1. The color factors are expressed in terms of 
the Casimir operators Cp = (A^ — l)/2Nc and Ca = Nc, where Nq denotes the number of 
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colors. The lengthy expressions for iVj and Ui in Eqs. (|^) and can be fonnd in |T^. For the 
gg process we have a slightly different way of splitting np the resnlts according to color than 
Ref. 0]. The choice in Eq. ensnres that the “Abelian” Uqed is identical to the QED part of 
19 QQ ill Rfif- [0 after taking into acconnt the usnal factor l/(2A"c') for replacing a photon 
by a gluon. Furthermore, compared to Ref. IQ an additional function Urf appears in Eq. (|^) 
since we are interested in the general case p/j ^ fip. Our unpolarized results fully agree with 
the corresponding expressions in 1.0 except for Uql in Eq. (6.22) of Ref. [^] which contains 
a numerically irrelevant misprint^. 

The real 2 —> 3 gluon bremsstrahlung corrections (R) can be split up in a similar way 
according to color structure and read 

1 




2(iV2 - 1) L 


{2CfYDqed + C^D 


A-^OQ 


+ D 


KQ 


( 6 ) 


and accordingly for the qq and gq processes |^. In order to isolate the divergencies of (Jg) 
appearing in the soft limit, which cancel the remaining singularities of the virtual matrix ele¬ 
ments, we “slice” the 2 —>• 3 result into a “soft” (S) and a “hard” (H) gluon part by introducing 
a small auxiliary quantity A 10. The kinematics of the soft gluon cross section is effectively 
that of a 2 —2 process, and all phase space integrations can be performed easily. Upon combi¬ 
nation with the virtual cross section the 1/e: IR and Xje^ IR-I-M singularities, all proportional 
to the n-dimensional Born cross section, drop out. Thus the “soft-|-virtual” (S-l-V) cross section 
becomes finite, except for 1/e: M singularities that cancel against “soft” x = 1 contributions, 
see Eq. (H), in the mass factorization procedure discussed below. Phase space integrations for 
the hard part are more subtle and require some care. As in we are interested in the double 
differential single inclusive cross section for the production of a heavy quark (or antiquark). 
For stable numerical simulations later on it is advantageous to perform the integrations over 
the phase space of the two not observed partons analytically as far as possible. To achieve 
this requires extensive partial fractioning to reduce all phase space integrals to a standard form 


1^, |T^, ^ . A sufficient set of four- and n-dimensional integrals are again conveniently collected 
in but we have recalculated and confirmed this set. Analytical results for the S-l-V cross 


sections will be given in ; the hnal expressions for the hard part are too lengthy but can be 
found in our computer code. 

To obtain a finite result, all remaining M singularities have to be removed by the standard 
mass factorization procedure which makes use of the fact that any collinear singular partonic 
cross section daf^ can be schematically written as 


daff{ki,k2)= / dxi / dx2Tii{xi, fiF/k‘)^mj{x2, ddim{xiki,X2k2) , 

Jo Jo 


( 7 ) 


^We thank J. Smith for his help in clarifying this issue. 
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where fci _2 are the momenta of the incoming partons. Here daim denotes the finite partonic 
cross section, /i is the scale introduced to render the coupling g dimensionless in n dimensions, 
and the transition functions Tij are given up to NLO by 


a.. 


Vifix,g.F/tj) = - x) + — 

ZTl 


5 ( 0 ) I 

ij 


X 


2 

-+lE 

e 


In dvr + In ^ ) +Cp (a;) 


( 8 ) 


with the LO splitting functions. Therefore all M singularities, Eq. (J^), can be absorbed uni¬ 
versally into the bare parton densities upon calculating physical, i.e., hadron-level, observables. 
The choice of the factorization scale yip and the arbitrary functions ^ij reflects the amount of 
arbitrariness in the separation of short-distance and long-distance physics. In the MS scheme, 
which we choose, one has ^ij = 0. 


Two subtleties, which show up in n-dimensionally regulated spin-dependent calculations be¬ 
yond the LO of QCD, have to be mentioned here. First and foremost, the projection operators 
onto certain helicity states, i.e., e^up^ and 75 , are of purely four dimensional nature, and there 
exists no straightforward and unique generalization to n 7 ^ 4 dimensions. We treat them by 
applying the internally consistent HVBM prescription where the Levi-Civita e-tensor con¬ 
tinues to be a genuinely four-dimensional object. 75 is defined as in four dimensions, implying 
75 } = 0 for /i = 0, 1, 2, 3 and [ 7 ^, 75 ] = 0 otherwise. The price to pay are {n — 4) dimen¬ 
sional scalar products, usually denoted by fc ■ p (“hat momenta”) appearing alongside the usual 
n-dimensional scalar products. In our case only a single hat momenta combination = —pPp 
appears in the polarized 2 —3 cross section and can be accounted for by an appropriately mod- 
ihed phase space formula |^. These contributions are inherently of 0(e) and only contribute 
to the final result when they pick up a 1/e pole. Secondly, the unphysical helicity violation at 
the qqg vertex in the HVBM scheme, which is reflected by Pqq^ 7 ^ APg® in n dimensions, has 
to be undone by an additional finite renormalization = — 4 C'i 7 ’(l —x) in Eq. (H) in the 

conventional MS scheme. Only then dAaqq = —d&qq is obtained, the result which is expected 
due to helicity conservation. 


3 Numerical Results and Phenomenological Aspects 


Before presenting results for hadronic heavy flavor distributions we first discuss the total par- 
tonic subprocess cross sections djj, i, j = q,q,g. They can be expressed in terms of LO and 
NLO functions and respectively, which depend only on a single scaling variable 

p = s/ (4m^) — 1: 










le-06 le-05 0.0001 0.001 0.01 0.1 1 10 100 1000 

ri=s/(4 m^) -1 

Figure 1: {vr? jc?^agg in NLO (MS) and LO as a function of// according to Eq. (1), where we have 
set jjLp = jjLji = m for simplicity and Anag = 2.7 as appropriate for charm production. 

where ag = as^jjip) and = (IIC'a —2n//)/3. Hence the are particularly suited for studying 
the main features of the NLO corrections in the most transparent way. The are non-trivial 
functions of rj and can be easily obtained from our double differential analytical results for the 
partonic cross sections by numerical integrations. In the unpolarized case they have been cast 
into a compact semi-analytical form |Q for fast numerical calculations of the total hadronic 
heavy flavor cross section which is a desirable future project in the polarized case. The can 
be derived just from mass factorization, since terms proportional to In originate only from 
Eq. (H). The last term in Eq. vanishes for the standard choice fip = fJ^R- In NLO this term 
follows straightforwardly from the LO result by replacing ag «« (^I + In thanks to 
the renormalization group invariance of the cross section. 

In Fig. I^we present the gluon-gluon subprocess cross section {w?/al)agg in LO and NLO for 
= hR = nr as a function of rj in the MS scheme. The threshold for QQ production, s = 4m^, 
is located at r; = 0. It turns out that the NLO corrections are signihcant in the entire rj range. 
At threshold the polarized and unpolarized cross sections are equal, thus Eq. implies that 
\Mgg \‘^{-\—) —> 0 as r] —0. Unlike in LO where agg approaches zero at threshold, it tends 
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to a constant in NLO, 2 {n^- i) ~ + f] “Conlomb singnlarity” 

present in the S+V part. It should be noted that in the threshold region logarithms from 
soft gluon emissions also contribute signihcantly even at the lowest t] shown. In the high 
energy domain, 77 —> oo, our polarized and unpolarized results behave rather differently. Here 
Feynman diagrams with a gluon exchange in the t-channel drive the unpolarized NLO result 
to a plateau value |Q whereas the polarized NLO cross section vanishes like the LO one, i.e., 
\Mgg\^{++) —>• \Mgg \^{-\—) iu Eq. 0.8 7] ^ oo. Similur observations have been made in the 
photoproduction case 'yg —> QQ |^. The scaling function fgg fulhls Afgg = —fqq after taking 
into account the additional hnite factorization mentioned above to restore helicity conservation 
at the qqg vertex. The behaviour of fqq was discussed in 0 and shall not be repeated here. The 
genuine NLO scaling function fgq is numerically much smaller than fgg and can be found in 
T 8 | but exhibits the same high-energy rj ^ oo behaviour as fgg, i.e., fgq approaches a plateau 


while ^fgq vanishes. 


The physical, i.e., experimentally observable, total cross section is obtained by convoluting 
the partonic cross sections in Eq. (P) with the appropriate flux of parton densities evolved to 
the scale yip, 

=y^l dxi dx2 fi{xi, fip) fj{x2, fir) , fip, , ( 10 ) 

' ^ I 4m2 / 4771^ 

ij S x-^S 

where S is the available hadron-hadron c.m.s. energy and s = X 1 X 2 S. In a similar fashion 
differential heavy (anti-)quark inclusive distributions like (Pd/dppdpL can be derived by con¬ 
volution with our double differential partonic cross sections. It also should be kept in mind 
that beyond the LO of QCD, parton densities and partonic cross sections have to be taken in 
the same factorization scheme in order to guarantee that Eq. (|T 0 |) is independent of unphysical 
theoretical conventions up to the order in considered in the calculation. 


One of the main motivations of performing the NLO calculations was to reduce the depen¬ 
dence on the choice oi pip and pLp which is completely arbitrary in LO and can lead to sizable 
ambiguities in predictions for d{S, m?) and the corresponding spin asymmetry to be dehned 
below. In Fig. ^we demonstrate that the NLO results for the polarized charm production cross 
section are indeed more robust under scale variations than LO estimates. In the left panel of 
Fig.0 we vary pp and pp independently of each other in the range p\ = Rw? and pf, = Fmf 
with 1 < R,F < 4.5 for hxed m = 1.4 GeV at a typical RHIC energy of '/S = 200 GeV using 
the GRSV “standard” set of polarized parton densities |^. In the right part of Fig. ^ we em¬ 
ploy the conventional choice pp = pp and vary pp and m = C GeV in a typical range for the 
charm pole mass. In order to better visualize the uncertainties due to scale and mass variations 
we show Aa{R, F,C)/Aa{R = 2.5, F = 2.5, C = 1.4) — 1, i.e., the deviation in percent of the 
total polarized charm production cross section according to Eq. (R) for variable pp^p and m 












Aa(R=HR/ m^, F=|Xp/ m^, C=m / GeV) / Act(2.5, 2.5, 1.4) - 1 



Figure 2: Deviation [in %] of the polarized total charm cross section in LO (dotted) and NLO 
(solid) from a reference choice ("O-pin" marker, see text) - left part: as a function of ftp and ftp 
for fixed m; right part: as a function of fxp and m with fxp = ftp, here the LO result is multiplied 
by a factor (-1). Corresponding contour lines in steps of 20% are given at the base of each plot. 

with respect to a reference cross section taken at fixed = 2.bm? with m = 1.4 GeV. 

To better guide the eye, contour lines in steps of 20% are plotted at the base of each plot. Here 
we also indicate the common choice /ir = ^p and m = 1.4 GeV (thin solid lines) in the left and 
right part, respectively. 

The NLO result in the left part of Fig. ^ is considerably “flatter” than the LO result 
with respect to variations of pi? but shows, however, slightly more variation with /ir. Not 
unexpectedly and more importantly it turns out that the usual choice fip = ^p almost coincides 
with the contour line of zero deviation from the reference cross section in NLO, in stark contrast 
to the situation at LO. This leads to the improved stability of the NLO prediction in the 
right panel of Fig. for variations of a common scale fip = fip at & given charm mass m. 
Here variations of the charm mass cause the major uncertainty of about ±30% in the NLO 
predictions. In LO we hnd considerable uncertainty stemming from variations of /if on top of 
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that. It should be also noted that qualitatively similar results are obtained for y/S = 500 GeV 
and bottom production at RHIC. Usually in NLO the terms proportional to and 

ln/i^//i|, in Eq. (|^) start to have a compensating effect for different choices for /ij? and /i/j and 
also provide some guidance that ///■ ~ 0 {m) and /ij? ~ /x/j in order to avoid large logarithms in 
the hard partonic cross sections. Ultimately one expects the dependence on /x^ and fiR to be 
reduced more and more if higher and higher orders in are considered. However, as was briefly 
explained above in connection with Fig. |^, in the reaction studied here new types of Feynman 
diagram topologies enter the calculation for the hrst time at the NLO level, whereas in next-to- 
NLO (NNLO) and beyond no qualitatively different diagrams appear. Hence in a sense NLO is 
the hrst “complete”, non-trivial order of perturbation theory for heavy havor production, and 
it is pleasing that scale stability improvements nevertheless clearly set in without considering 
NNLO corrections which seem unattainable at this time. 


Instead of measuring polarized cross sections like Aa{S,m‘^) directly, experiments will usu¬ 
ally study the related longitudinal spin asymmetry dehned by 

Aa{S, m?) 




( 11 ) 


(j(S', m?) 

in case of the total cross section and accordingly for differential heavy quark distributions. The 
experimental advantage of this quantity is that one does not need to determine the absolute 
normalization of the cross section Aa{S,m‘^) which is usually difficult to obtain. However, one 
should keep in mind that the situation in the unpolarized case is far from clear, in particu¬ 
lar concerning bottom, and hence it would be reasonable to determine the unpolarized and 
polarized cross section separately. We note that for small variations of the scales the relative 
deviation of the asymmetry can be written as ^ ~ H turns out for the variations 

of hf, Hr, and m considered above that and ^ are almost equal in NLO, whereas they 
can differ strongly in LO. As a result it is even more true for the asymmetry that NLO results 
are highly stable, whereas the LO uncertainty is huge, in particular for the choice fiR = hF- 
We will explore this in detail in |^, but wish to point out here that LO determinations of Ag 
using the asymmetry alone will necessarily have a prohibitively large theoretical error, a NLO 
analysis is a must in that case. 


Finally, let us turn to the important question of whether heavy flavor production at RHIC 
can be used to discriminate between different polarized gluon densities. To address this ques¬ 
tion thoroughly one has to take into account an estimate of the statistical signihcance of a 
measurement of a heavy quark spin asymmetry at RHIC. Compared to direct photons or jets 
which are directly observed in the detector this is a rather involved problem for heavy flavors. 
With the PHENIX detector at RHIC charm and bottom quarks can be identihed only through 
their decay products, preferably leptons. However, the electron and muon detection is rather 
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limited in pseudo-rapidity, \r]f>\ < 0.35 and 1.2 < \r]^\ < 2.4, respectively, and cuts in the lepton 
Pt have to be imposed in order to separate charm and bottom. Since heavy flavor decays usually 
have a multi-body kinematics and may proceed through cascades, cuts on the observed leptons 
are difficult to translate back to the calculated parton, i.e., heavy quark, level. One possibility 
is to rely on Monte Carlo simulations of heavy quark decays, for instance, on PYTHIA 


which are quite successful and tuned to a wealth of data. PYTHIA can be used to generate 
“efficiencies” for observing a heavy quark within a certain bin in pt and p in the detectors at 
RHIC. If properly normalized to the total number of heavy quarks generated in that particular 
bin, £eff should become independent of all the details of the heavy quark production mechanism 
assumed in PYTHIA. This is essential for an unbiased extraction of Ag. 

Exploiting this idea, a first numerical study for the PHENIX detector has been performed^ 
The resulting efficiency ees{pT,v) for a charm quark produced with transverse momentum pt 
and pseudo-rapidity p to be detected via its decay electron anywhere in the PHENIX acceptance, 
with the electron trigger allowing p^ > 1 GeV, is approximately given by 

-9.79 + 4.58(pT/GeV)i'^ 


ees{pT,P]PT > 1 GeV) = C exp 


with 


(pT/GeV)i-73 + 1.74 

C = exp \p/ (4.06 exp [— (pr/1.05 GeV)°''^^]) (PT/2.48GeV) 


( 12 ) 


A prediction for the charm cross section as measurable with PHENIX is then obtained by 
convoluting our double differential partonic results with SeS in Eq. ([T^), 

(Pa 




o-eff(PT > 1 GeV) = / dpT / dr/ £efr(PT,^?;PT > 1 GeV) , 

Jo J -ri^ax dpTdfj 


(13) 


where = l\/S — AnP and p^^^ = — | In ^ appropriate kinematical 

2'^ I 2 i+^i_4p2/(5_4m2) nn n 

limits. In the future |]^ we plan to use a set of efficiencies with several different cuts on, or 
bins in, p?p to generate different des- For the time being, different cuts in p^ are simulated by 


limiting the charm transverse momentum px instead, while still using Eq. (pY]), i.e., 
deff(PT > Pt”) - 

This expression has been used for the results shown in Fig. 


dpT I dp EesipT, P]Pt> ^ GeV) 


(14) 


To study the sensitivity of the charm production asymmetry at RHIC to Ag, we use a range 
of recent [^, |^, ^ and old |^, ^ helicity densities in Fig. |^. These sets mainly differ in the 
assumptions about Ag. Note that for calculating the required unpolarized a in A = Aa/a we 
have used in each case the underlying set of helicity averaged parton distributions as specified 
^We are grateful to M. Grosse Perdekamp for providing these efficiencies. 
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Figure 3: The NLO charm asymmetry A at a/^ = 200 GeV for PHENIX at RHIC as a function 
of using Eq. (|5). For a better separation of the curves A is rescaled by 

Recent and old sets of helicity densities are distinguished by thick and thin lines, respectively. Also 
shown is an estimate for the statistical error using a luminosity of £ = 320 pb“^ (see text). 


m 


24, 


i^, For consistency, m is also taken as in these fits, i.e., m = 1.4 GeV for the 
modern and m = 1.5 GeV for the old spin-dependent densities. All results are obtained for the 
choice /ip = /i|. = 2.5(m^ +Pt)- If is immediately apparent from Fig. that charm production 
at RHIG can be very useful in pinning down Ag. The estimated statistical error for such a 
measurement, 6 A = ^- 7 ==, assuming a luminosity of £ = 320 pb“^ and a beam polarization 
of Pp ~ 0.7 [0 is signihcantly smaller than the total spread of the predictions. This is true in 
particular if we take the large GS A Ag into account, which is still compatible with recent 
data. The GS A asymmetry is so large that it had to be scaled down by 0.7 to £t well into the 
same plot. Note that very small gluons, e.g., the oscillating Ag of GS G |2^, are at the edge 
of being detectable. There are three groups of helicity densities which will be indistinguishable 
within the errors by charm production: GRSV’OO std. & DS i-, AAG 1 & 2 , and DS iii-|- & 
GRSV’96 std. & val. Here the gluon densities are too similar in the range of x predominantly 
probed by this process. We will map the range in x where Ag is accessible by heavy flavor 
production at RHIG in detail in 
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Figure 4: Ratio of the asymmetries in NLO and LO, 74 (NLO)/y4(LO), with 74 (NLO) as shown in 
Fig. 1^. The ratio of unpolarized cross sections (j(LO)/cr(NLO) used in the calculation of A is also 
shown for comparison (thin solid line). 


Finally we take a look at the importance of the NLO corrections for the heavy flavor spin 
asymmetry. Fig. ^ shows the ratio of the charm spin asymmetries calculated in NLO (as shown 
in Fig. ^ and LO for the range of modern helicity densities [^, ^ 


One can infer that 


the NLO asymmetries are generally smaller than the LO ones by a factor of about three. In 


the case of the AAC helicity densities we find an even larger suppression. Furthermore, 
there is considerable dependence on which will inhibit the use of constant “K-factors” 
to estimate the NLO results from LO ones. It should also be pointed out that much of the 
reduction of the asymmetry in NLO stems from the unpolarized cross section in NLO, which 
is about a factor two larger than the corresponding LO result. This is shown by the thin solid 
line in Fig. ^representing the ratio ct(LO)/(t(NLO) obtained with the GRV’98 parton densities 


311. The sizable difference of the asymmetry predictions in LO and NLO means that the LO 


and NLO gluon helicity densities extracted from a future asymmetry measurement will differ 
considerably. Whether this will be consistent with data from other processes has to studied in 


a global analysis, e.g., along the lines suggested in Ref. |30 
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Further studies of the uncertainties and predictions for bottom production will become avail¬ 


able soon [Q. There we will also present more details concerning the calculational techniques 
that have been used and analytical results for the matrix elements that we have obtained. 
In case one wants to extend onr glnon-glnon resnlts to the prodnction of glninos, one has to 
take into acconnt the following colorfactor replacements in Eqs. (^) and (|^), which already 
inclnde the color-average: replace the prefactor , I . —> and set Cp = Ca = Nc inside 

the sqnare brackets. Before doing the latter, one has to nse the identity 1 = C\ — 2CfCa 
for the KQ parts, e.g., Urq = {C\ — 2CfCa)Ukq —*■ —NqUkq- We note that one conld 
nse this identity anyway to rewrite the resnlts with the same nnmber of colorfactors, e.g., 
C\UoQ + Ukq = C\{UoQ + Ukq) — 2CfCaUkq, bnt this has nnmerical disadvantages dne to 
the length of the expressions for Uqq and Ukq- For a complete NLO snpersymmetric QCD 
calcnlation of spin-dependent glnino prodnction one wonld have to take into acconnt that the 
rnnning of the strong conpling will be changed and that other snbprocesses will contribnte. 
However, as always, glnon-glnon scattering is expected to be the dominant snbprocess. 


4 Summary 

To snmmarize, we have presented the hrst complete NLO QCD calcnlation for the spin- 
dependent hadroprodnction of heavy qnarks. The NLO resnlts have considerably less nncer- 
tainties stemming from variations of the nnphysical factorization and renormalization scales. 
They even become nearly independent of the scales when the conventional choice /Xi? = is 
employed. There remains, however, an nncertainty of Aa of abont ±30% from variations of the 
charm mass within a reasonable range. We have made predictions for the charm asymmetry 
that can be measnred soon with the PHENIX experiment at RHIC. These predictions inclnde 
an efficiency fnnction which describes hadronization, decays of the heavy qnarks, experimen¬ 
tal cnts, and the detector geometry. In this way it was demonstrated clearly that Ag can be 
constrained considerably with heavy flavor prodnction, even if cnrrent experimental limitations 
are taken into acconnt realistically. Finally, we have shown that as in the nnpolarized case, 
LO calcnlations cannot be snbstitnted in any simple manner for the fnll NLO resnlt. Onr 
work is of particnlar relevance also to the “heavy qnark enigma” arising from discrepancies of 
nnpolarized bottom, bnt not charm, prodnction data with theory. On the one hand one now 
will be able to compare spin-dependent KHIC data for charm and bottom prodnction to a NLO 
prediction. On the other hand onr calcnlation provides the major part of the spin-dependence 
of a possible “new physics” explanation of this discrepancy in terms of snpersymmetry by a 
simple adjustment of colorfactors. More details and fnrther phenomenological stndies can be 
fonnd in Q. 
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